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Abstract 

We apply boson expansion methods to an extended Lipkin-Meshkov-Glick 
model including anharmonicities in analogy with previous microscopic calcu- 
lations. We study the effects of different approximations present in these cal- 
culations, among which the truncation of the hamiltonian and of the space, in 
connection with the study of the properties of two-phonon and three-phonon 
states. By comparing the approximate results on the spectrum with the exact 
ones we conclude that the approximations made in the microscopic calculations 
on two-phonon states are well justified. We find also that a good agreement 
with the exact results for the three-phonon state is obtained by using a bosonic 
hamiltonian truncated at the fourth order. This result makes us confident 
that such approximation can be used in realistic calculations, thus allowing a 
theoretical study of triple excitations of giant resonances. 
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1 Introduction 



Collective excitations have been known for many years in nuclear physics 0, both in 
the low-lying spectra and in the Giant Resonance (GR) region. A basic microscopic 
theory for such collective modes is the Random Phase Approximation (RPA) 0, |3| 
which can be seen as the lowest order in a boson expansion such that the hamiltonian 
can be put in the form of a sum of hamiltonians of harmonic oscillators, each one 
corresponding to a collective mode (phonon). Therefore, RPA predicts the existence 
of one-phonon, two-phonon,...etc states with a harmonic spectrum. In addition to the 
well known low-lying two-phonon states, recently heavy ion inelastic scattering exper- 
iments at intermediate and relativistic energies and double charge-exchange reactions 
have shown the existence of states in the high excitation energy region which can be 
described as a GR built on top of another GR j|, |5| . The study of the properties of 
these states allows to test our comprehension of the GR's as small amplitude vibra- 
tions and therefore the harmonic picture. The systematics on the energies and widths 
is in qualitative agreement with the harmonic approximation, namely the energy of a 
doubly excited GR is nearly twice that of the single GR and its width is between \pl 
and 2 times larger. However, the inelastic cross sections, when calculated within the 
harmonic picture, are almost always smaller than the measured ones. In particular, 
the data concerning Coulomb excitation show a discrepancy ranging from 30% up 
to a factor 4, according to the nucleus studied. In order to understand the origin of 
this discrepancy, corrections to the harmonic approximation have been proposed || 
by including anharmonicities in the internal hamiltonian and non-linearities in the 
external field. As shown in |J, small anharmonicities in the excitation spectrum of 
the target nucleus can lead to a large enhancement of the Coulomb excitation cross 
section. The model used there was an oversimplified one, namely the target was 
described as an anharmonic linear oscillator. The parameters of the cubic and the 
quartic terms in its hamiltonian were fixed so that the energy of the second excited 
state was shifted down by ~ 2 MeV with respect to twice the energy of the first 
excited state. In [|7| a 3-dimensional extension of this model was considered and sim- 
ilar effects coming from anharmonicities were found. We remark that, in the case 
of Coulomb excitation of 136 Xe in the reaction 136 Xe + 208 Pb at E/A = 700 MeV, 
the experimentally observed peak, interpreted as the Double Giant Dipole Resonance 
(DGDR), is shifted by ~ —2 MeV from that expected in the harmonic case ||. A 
more refined study was performed in where anharmonicities and non-linearities 
were included by starting from RPA and extending it by means of boson mapping 
techniques ||, Such microscopic approach was applied in || to realistic cases. In 
the case of 208 Pb + 208 Pb at E/A = 641 MeV, for example, it was found that the 
inclusion of anharmonicities and non-linearities gives rise to an enhancement of 30% 
of the cross section in the region of the DGDR, bringing the theoretical results closer 
to the experimental ones. Important contributions coming from other two-phonon 
states in the same energy region were also found. 

The studies reported in the above quoted papers have of course some limitations. 
The model used in |J is clearly very schematic and has, among others, the drawback 
that it is based on a purely bosonic description so that Pauli blocking effects are not 
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included. The approach presented in || is certainly much more realistic and the effects 
of the Pauli exclusion principle are taken into account to some extent. However, for 
computational reasons, only one- and two-phonon states were considered there. 
The purpose of this paper is to present an analysis based on an extension of the 
Lipkin-Meshkov-Glick (LMG) model [D]] including a residual interaction between the 
phonons, in analogy with the microscopic approach used in ref.Q. This model has 
several advantages with respect to the others. Because of its group structure, this 
model is exactly solvable. Besides, it directly takes into account the Pauli principle. In 
this context, we discuss several approximations, corresponding to the cases considered 
in H 3 and, by direct comparison with the exact results, we test how severe are the 
limitations of these approximations. We focus on one hand on two-phonons states, 
studied in ref.0 and on the other hand on the properties of three-phonon states. In 
fact, in order to have a more stringent test on the validity of the harmonic picture, 
experimental and theoretical studies of triple excitation of GR's should be envisaged. 
From the theoretical point of view, the same approach as the one in can be used, 
but this would imply huge calculations. In order to make them feasible in the context 
of boson expansion methods, one may consider an approximation in which the same 
fourth order hamiltonian used in |J is diagonalized in a space containing up to 3- 
phonon states. We present a test on how well this approximation works, in the context 
of the extended LMG hamiltonian. We will see that the exact results are very close 
to the approximate ones in all the cases studied for both the second and the third 
excited states. 

The paper is organized as follows. In section 2 the LMG model is shortly reviewed 
and our extension introduced. In section 3 the boson mapping of this hamiltonian 
is presented. In section 4 the approximate results on the energy of the first, second 
and third excited states, obtained by truncating at different levels the mapping, are 
compared among themselves and with the exact calculations. 



2 The model 

2.1 The Lipkin-Meshkov-Glick model 



In the original Lipkin-Meshkov-Glick (LMG) model |TTJ a finite number Q of particles 
can be arranged in two levels separated by e in energy. Q different quantum states are 
available in each level. Each particle is therefore identified by two quantum numbers. 
The first a indicates if the particle is in the upper [o = +) or in the lower (a = —) 
energy level. The second quantum number s is related to the Q different available 
quantum states in each level. We will consider systems whose hamiltonian can be 
written as function of the operators K + , K_ , K Q : 
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5=1 

= -XK a +, s a +,s- a -, s a -J 

1 8=1 

where aj. s (and the hermitian conjugate a CT)S ) creates (annihilates) a particle in the 
quantum state (a, s). Since the operators K + , K_ , K satisfy the SU(2) commutation 
relations : 

[K+, KJ\ = 2K [K , K ± ] = ±K± (2) 

they are often called the quasi-spin operators. According to (^]), there can not be two 
particles in the same quantum state, that is the Pauli exclusion principle is included 
in the model. A system of O particles has 2 Q states. However, since the hamiltonian 
is a function of the generators of SU(2) and K 2 is a Casimir operator of such algebra, 
the space can be separated into subspaces, each corresponding to an eigenvalue K. 
The state with all particles in the lower level, a = — , is an eigenvector of K 2 and Kq, 
belonging to their maximum and minimum eigenvalue, respectively. This state, with 
the ones generated by applying K + to it, are the elements of a subspace of dimension 
O + 1, that we will denote by |0/2,m) with m G [—0/2,0/2]. 
Following ([!]) the original hamiltonian of LMG [[I I]] can be written as : 

H LMG = eK + V X K+K_ + V 2 (K + K + + K_K_) (3) 

where 

K = K + 0/2 (4) 

The eKq term can be viewed as corresponding to the Hartree-Fock part. The lowest 
eigenstate of K corresponds then to the uncorrelated HF ground state |0/2, —0/2). 
In order to stress the analogy with the microscopic calculations in [Q, we will often 
denote by h (hole) a (— ,s) state, i.e. a single particle state which is occupied in 
|0/2, —0/2), and by p (particle) a (+, s) single particle state. Therefore, by looking 
at eq.(|lD, we can say that the residual interaction in (^) only contains particle- hole 
terms a) p a} p ,ahCLh' and a^a^ap/ah' terms, that is those usually included in RPA. 
If Vi = V2 = then Hlmg = sKq. The corresponding eigenvalues are ne with 
n G [0, 0] . From now on we will denote these states by \n >. The energy eigenvalues 
are equidistant and form a harmonic spectrum truncated at + 1 levels. If V\ or 
V2 are different from zero, then the energy spectrum is not harmonic any more and 
only if V2 7^ the energy eigenvectors correspond to a mixing of |n) states. In fact, 
V\K + K_ is diagonal in the \n) space. This term of ph — p'h' type does not mix 
states with different particle-hole numbers. On the contrary, the pp' — hh' type term, 
V2 (K+ + K 2 ^ mixes states with n and n±2 particle-holes. To conclude this part, we 
would like to recall that the LMG model is exactly solvable using group techniques 
and that it includes the Pauli exclusion principle. 
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2.2 An extension of the LMG model 



The LMG model in its original form already includes some anharmonicities, essentially 
those related to the fact that the Pauli principle is treated exactly. However, parts 
of the residual interaction are neglected, namely the pp' — p"p"', pp' — p"h, hh' — h"p, 
hh' — h"h"' ones, which were considered in the microscopic calculations |J where the 
anharmonicities arising from them were also studied. In order to simulate them, we 
propose an extension of the LMG hamiltonian which is still quadratic in the K + , K_ 
and K operators : 

H = Hlmg + AV = H lmg + V 3 (K+K + KqKJ) + V 4 (K - 1)K (5) 

The K + Kq term and its hermitian conjugate introduce a coupling between \n) and 
\n ± 1) whereas the last term shifts the energies of the |n) states, except those with 
n = and 1. Therefore, the eigenstates \<f> a > of the hamiltonian (|]) are superposi- 
tions of the \n > states 

|0 a >=]>»> (6) 

n 

In section 4 we will compare the exact eigenvalues of the hamiltonian (||]) with 
those corresponding to the bosonic hamiltonian obtained by mapping the fermionic 
one up to the fourth order. 



3 The boson hamiltonian 

We apply boson expansion methods to the fermionic hamiltonian (|5]) and truncated 
the so obtained boson hamiltonian to second and fourth order. In ref.sfl^, 13 1 a 



similar study was done for the original LMG hamiltonian (Q). In the present case, the 
quadratic hamiltonian corresponds to RPA while the quartic can be directly compared 
with what was done in 0, where some further approximations were performed, as 
discussed below. 

Let's take the normal ordered Holstein-Primakoff boson expansion of the SU(2) gen- 



erators 12, 14, 15 



(Jf-) 6 = E£ fli(6 t )'6 m (7) 
(Ko)b = tfb 

where 



E 7 1 ir (n 



m 



1/2 



In the case of SU(2), the Hage-Hassan-Lambert mapping used in ref.0 is equivalent 
to the Holstein-Primakoff mapping (|7|,[8D |16|] with the shifted operator (Kq)^ . By 
this shift one eliminates terms linear in W and b in the hamiltonian which corresponds 
to a redefinition of the mean field. 

Let us then consider the bosonic mapping of our extended hamiltonian ([|) in the 
particular case V2 = V\/2 which is the value used in the numerical applications 
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we present in the next section. By including all the terms of the expansion (^) 
which contribute up to the second order in the tf, b operators, we get the quadratic 
hamiltonian 

H?> = e(l + 6 1 pb + 6 f(l~ ^) V2 ^ + b 2 ) (9) 

where 

5 1 = QVi/e (10) 
is related to the strength of the particle-hole residual interaction. The hamiltonian 

(2,) 

HI can be written in diagonal form by introducing the Bogoliubov transformation 

&t = xgt + Y Q 
b = XQ + ygt 

and imposing 

#20 = (12) 

with the condition 

X 2 - Y 2 = 1 (13) 
which guarantees the correct commutation relation 

Q,Qt] = i ( M ) 

Thus one gets the RPA hamiltonian 

= HuQ ] Q (15) 



1 \ 1/2 

H n = e{l + 5i)(X 2 + Y 2 ) + 2XYS 1 e 1 - - (16) 



where 

= P n -+- ^,vy 2 + y 2 ^ f i - 

The conditions ( II^LO ) determine the X and F's amplitudes as solutions of the set of 
equations : 

X 2 -F 2 = 1 

e(l + 5 1 )XF + ^(l-i) 1/2 (X 2 + r 2 ) = (17) 



Of course the terms in V 3 and V4 do not appear at the quadratic order. In eq.(15), 
the H 00 which corresponds to a shift of the eigenenergies was omitted. 
The violations to the Pauli exclusion principle introduced by truncating the boson 
expansion at the lowest order as well as a contribution coming from the residual 
interaction of the pp' — p"p"', pp' — p"h, hh' — h"p, hh' — h"h!" type can be partly 
introduced by going one step further, i.e. including all terms up to the fourth order, 
as it was done in ref.||. At this level, the terms in V3 and V4 will therefore enter. 
Their presence, together with the corrections for the Pauli principle, will give rise to 
an anharmonic bosonic hamiltonian. 

By mapping the hamiltonian of eq.(B) up to fourth order we get: 
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H?> = e(l + 8 x )Wb + 5 -f(l- ^) 1/2 (fet 2 + 6 2 ) 

+v 3 VU (1 - 1) 1/2 (&+ 2 6 + b% 2 ) + (v 4 - bH 2 

(bH + tfb 3 ) 



+ 5 -f{l 



1/2 



1/2 



which we rewrite in terms of the and Q operators as : 



Hl A) = FnQtg + # 30 (Qt3 + q 3) + # 21 (gf2g + gtg 2)+ 
H 31 (Q^Q + QtQ 3 ) + H 22 QVQ 2 + H m {Q^ + g 4 ) 



(19) 



whose coefficients are given in the appendix (eqs.^T]-[27]), together with the new equa- 
tions for X and Y (|20|) coming again from the conditions ( |il?i|l3l) . In eq. (|T9|) the 
Hqo was omitted as well as a term, linear in the and Q operators, which would 
introduce a redefinition of the mean field, in analogy with what was done in |J. 
The quartic hamiltonian (|19|) corresponds to that used in |J where, however, some 
approximations were introduced in order to make feasible the calculations in the real- 
istic cases considered there. First, only one- and two-phonon states were considered. 
Therefore, the terms in H 30 , H 3 i and were not effective. Second, the X and Y 
amplitudes appearing in the fourth order hamiltonian were not recalculated but taken 
equal to those obtained at the second order, i.e. the RPA ones. In order to get an 
indication on how good these approximations on the space, on the hamiltonian and 
on the X and Y amplitudes are, in the next section we will study them within the 
present schematic model. 



4 Results and discussion 

First of all we have to fix the parameters entering in the hamiltonian @. For the single 
particle energy we use the parametrization e = Al/A 1 ^ 3 MeV. We take V 2 — 0.514 
and the strength V\ is fixed so that the first excited state lies at an energy around 
80/y4 1//3 MeV corresponding to the systematics of GDR in nuclei. This criterion gives 
V\ = 1.2 MeV. We have studied the behaviour of the energies of the three lowest 
states, Ei, E 2 and E%, as a function of V4 for two values of V3, namely V3 = MeV 
(Fig.l) and V 3 = 0.25 MeV (Fig.2). The latter value for V 3 gives < 2|AV|1 >^ 1 MeV 
in analogy with the microscopic calculations [H. Note that the sign of V3 is irrelevant. 
As far as the sign of V4 is concerned, we show results only for negative values, which 
give a downward shift of E 2 and E 3 with respect to 2Ei and 3Ei, respectively, i.e. 
the harmonic (RPA) values. In the figures we show the results obtained for Q = 8. 

To compare the spectrum of the exact hamiltonian and its different boson expansions 
we study the energy differences (E n — E n ^i) between the lower excited states. In RPA 
they have the common value 15.11 MeV for the considered parameters, independent 
of n, V3 and V4. In the figures this value is represented by a dashed line. This should 
be compared with the exact results for the fermionic hamiltonian eq.(5), shown as 
solid lines. We see that, even for V 3 = V4 = 0, the RPA results deviate from the exact 
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-0.4 -0.2 0.0 

V 4 (MeV) 

Figure 1: Energy differences as functions of V4 for a fixed value of V3 = MeV. 
The ones corresponding to the fermionic hamiltonian (g) are plotted as solid lines. 
The others correspond to different approximations on the bosonic expansion of the 
hamiltonian (see text). 

ones. Namely, the first three differences obtained in the exact calculations are equal 
to 15.62 MeV, 14.38 MeV and 13.03 MeV, respectively. 

The agreement with the exact energies can be improved by using the quartic hamilto- 
nian, eq.(19). When the bosonic hamiltonian is diagonalized in the space containing 
up to two-phonon states we get the results presented in Fig. 1 and 2 by dotted lines. 
We remark that now the first excited state is very close to the exact one, while for 
the second one there is a discrepancy of not more than 250 KeV. In the more realis- 
tic microscopic calculations presented in 0] a further approximation was introduced. 
Instead of reobtaining the X and Y amplitudes by solving the equations analogous 
to eqs (20), the RPA amplitudes were used. In the present model, this corresponds 
to use in eq. (19) the solutions of eqs (17) rather than those of (20). The results 
obtained are indistinguishable from the dotted lines in Fig. 1 and 2. This supports 
the validity of the procedure used in ||. 
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-0.4 -0.2 0.0 

V 4 (MeV) 

Figure 2: As in fig. 1, but for V 3 = 0.25 MeV. 

In order to test how much the truncation of the space affects the results, let us now 
enlarge the space up to three-phonon states. Of course, in this space, the quartic term 
if 40 in (19) does not play any role. In the figures, we show the results corresponding to 
the complete quartic hamiltonian of eq . (|i~9"D , with the X and Y amplitudes solutions 
of eq.s (20) (dot-dashed lines). Comparing the energies of the first two states we 
see that the agreement with the exact ones is now almost perfect. In the enlarged 
space we can also study the third excited state which, in the harmonic limit, would 
correspond to a three-phonon state. In this case the approximate results differ from 
the exact ones by 1 MeV, which can be considered as a good approximation since 
the excitation energy of this state is about 40 MeV. We would like to point out 
that the results obtained in the same model space but corresponding to the quartic 
hamiltonian of eq (19) without the H 3i and H 3Q term and with the use of the RPA 
X and Y amplitudes, solution of eqs (17), (represented in figure 1 and 2 like solid 
circles) are practically coincident with the dot-dashed line. The use of an analogous 
approximation in realistic microscopic calculations would make them much simpler. 
In summary, we have seen that quite good results can be obtained for both the second 
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and the third excited state by diagonalizing the fourth order bosonic hamiltonian with 
coefficients calculated with the RPA A and Y amplitudes, in the space containing 
up to three-phonon states. Nowadays there is interest in searching for triple giant 
resonance excitations [T7j. The present results indicate that a realistic theoretical 



study of such states will be feasible in the near future. 



5 Conclusions 

We have applied boson expansion methods to an extended Lipkin-Meshkov-Glick 
model including a residual interaction between the excited states in analogy with 
the anharmonicities introduced in the microscopic calculations of ref.||. We have 
studied the effect of truncations of the space and of the hamiltonian as well as other 
approximations present in those calculations by comparing the approximate results 
on the energies of the first three excites states among them and with the exact ones. 
From the analysis presented we can conclude that the approximations made in || 
are well justified if one wants to study one- and two-phonon excitations. What about 
three-phonon states? Of course, in the context of boson expansion methods, the study 
of these states would need a boson expansion up to sixth order, but this would be 
a formidable task. The results shown make us confident that a quartic hamiltonian 
as that used in ||, diagonalized in an enlarged space including up to three-phonon 
states, may give reasonable results also in realistic cases. 
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Departamento de FAMN of the University of Sevilla where part of the work has been 
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supported by the Spanish DGICyT under contract PB95-0533-A, by the Spanish- 
Italian agreement between the DGICyT and the INFN and by the Spanish-French 
agreement between the DGICyT and the IN2P3. 



6 Appendix 



The A and Y coefficients of the Bogoliubov transformation obtained by imposing 
conditions (|l2]Jl6|) for the quartic expansion of the hamiltonian (|5]) are solutions 
of the equations 



X 2 -Y 2 

e(l + 6 X )XY + 5 -f (l - i) V2 (A 2 + Y 2 )+ 
[ {V 4 - Vi) XY (A 2 + 5F 2 ) + 3A (3A 2 F 2 + Y A ) 



(20) 







The coefficients of the boson hamiltonian are given by the following expres- 
sions 
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H u = e(l + 5i) (X 2 + F 2 ) + 2XF5xe (l - ^) 1/2 + 
4 (2X 2 F 2 + F 4 ) (V4 - V x ) + 6XYX (X 2 + 3Y 2 ) 

H 30 = v 3 xy (x + Y)Vnn- 

#21 = Ksv 7 ^ 1 1 — — J (x 3 + 2XF 2 + 2X 2 F + Y' 



and 



H. 



40 



K + 1/ 4 ) X 2 F 2 + AXF (X 2 + Y 



H 31 = (-V 1 + V 4 ) 2XY (x 2 + Y 2 ) + A (x 4 + 6X 2 F 2 + r 
H 22 = {-Vx + V4) (X 4 + 4X 2 F 2 + F 4 ) + QXY A (x 2 + Y' 



In particular, if Vi = V3 = V4 = 0, the relations 
reference [[L2[ . 



(2 

coincide with those 
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